Linear and nonlinear interactions in the dark sector 
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We investigate models of interacting dark matter and dark energy for the universe in a spatially flat 
Friedmann- Robertson- Walker (FRW) space-time. We find the "source equation" for the total energy 
density and determine the energy density of each dark component. We introduce an effective one- 
fluid description to evidence that interacting and unified models are related with each other, analyze 
the effective model and obtain the attractor solutions. We study linear and nonlinear interactions, 
the former comprises a linear combination of the dark matter and dark energy densities, their first 
derivatives, the total energy density, its first and second derivatives and a function of the scale 
factor. The latter is a possible generalization of the linear interaction consisting of an aggregate of 
the above linear combination and a significant nonlinear term built with a rational function of the 
dark matter and dark energy densities homogeneous of degree one. We solve the evolution equations 
of the dark components for both interactions and examine exhaustively several examples. There 
exist cases where the effective one-fluid description produces different alternatives to the ACDM 
model and cases where the problem of coincidence is alleviated. In addition, we find that some 
nonlinear interactions yield an effective one-fluid model with a Chaplygin gas equation of state, 
whereas others generate cosmological models with de Sitter and power-law expansions. We show 
that a generic nonlinear interaction induces an effective equation of state which depends on the scale 
factor in the same way that the variable modified Chaplygin gas model, giving rise to the "relaxed 
Chaplygin gas model" . 
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I. INTRODUCTION 



Certain high precision astronomical observations sug- 
gest that the Universe entered an accelerated expansion 
stage when the value of its scale factor was approximately 
a half of the current one. This important discovery, which 
was based on the observations of the brightness of a class 
of supernovas (SNIa) [l[ , has been confirmed by precise 
measurements of the spectrum of the cosmic microwave 
background (CMB) anisotropics 0] as well as the baryon 
acoustic oscillations (BAO) in the Sloan digital sky sur- 
vey (SDSS) luminous galaxy sample 0. This discovery 
has transformed Cosmology into a very active area in cur- 
rent Physics and will surely fix the bases of important 
advances in the future because the consensus between 
cosmologist points in the direction that the understand- 
ing of the phenomenon will probably require a unified 
comprehension of the gravitational and the other funda- 
mental interactions. The above aspect of the expansion 
of the Universe, which becomes manifest on very large 
scales, will be only detectable from very distant astro- 
nomical objects. 

Basically, the evidences indicate that the Universe was 
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dominated by nearly pressureless dark matter in the long 
lasting initial stage that goes back almost to its most 
early stages after the Big Bang. This was an epoch char- 
acterized by an every time slower expansion. However, 
the behavior got reversed and the Universe began an ac- 
celerate expansion under the domination of its dark en- 
ergy component characterized by a negative pressure [IJ . 
This behavior has lasted till present and most likely con- 
tinue for ever. At large scales, there are strong evidences 
for a spatially flat and accelerating universe transiting 
from a scenario dominated by matter accumulated by 
purely attractive gravitational effects to another domi- 
nated by a dark component dispersed by repellent gravi- 
tational effects. The reasons, causes and details of when 
this transition happened, it was not still understood. 
These are part of the queries to answer in those projects 
which are framed in this active investigation area inside 
the Cosmology. It is interesting to note that for a suffi- 
ciently intense acceleration one can speak rather of super- 
acceleration, then the possibility exists that the universe 
has a catastrophic end with a suddenly future singularity 
at a finite time (Big Rip) and a total disintegration of the 
well-known structures. 

To investigate the mechanisms that govern the dynam- 
ics of the evolution of the universe from its early stage 
until its recent accelerated phase we will consider fun- 
damentally two types of model denominated respectively 
interacting and unified models. The aspects of the evo- 
lution in which we will be concerned arc those connected 
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to theoretical descriptions of the dark matter and dark 
energy. 

In the interacting models, the source of Einstein equa- 
tions which describe the dynamics of the universe at large 
scale includes an aggregate of different material fluids 
and scalar fields that are conserved individually or in- 
teract among them. This is in principle the simplest, 
and perhaps the most obvious hypothesis, and it is in 
fact, the one that has provided more advances in the 
knowledge of thephenomenon of the recent acceleration 
of the universe |f|, @, and references therein. Follow- 
ing observational evidences we will consider three fun- 
damental components: baryons, dark matter and dark 
energy. Given the dynamical similarity between baryons 
and dark matter, we will make a simplified model re- 
placing both components with a nearly pressureless dust, 
while the dark energy will be described by a fluid with a 
linear equation of state. This will allow us to focus our 
investigations on models of two fluids with energy trans- 
fer. Our goals will be the following ones: on one hand 
we will mainly investigate linear and nonlinear interac- 
tions, they will be considered as functions of dark matter 
and dark energy densities, their first derivatives, the to- 
tal energy density with its derivatives up to second order 
and the scale factor. On the other hand we will analyze 
the relation between interacting and unified models. We 
will also examine the problem of coincidence: are the 
proportions of matter accumulated by gravitative effects 
and dark energy comparable at the present time for a 
strange coincidence or for a fundamental reason? Many 
models in the literature Q have been proposed to allevi- 
ate the problem of coincidence, for instance quintessence, 
k-essence, phantom, quintom, tachyon, etc. 

In unified models, the Einstein equations will have a 
single component working as dark matter and dark en- 
ergy at different stages. It interpolates smoothly between 
a matter dominated phase in the early stage and dark 
energy in the late stage of the evolution, so inducing an 
accelerated expansion of the universe. Consistently, the 
universe evolves from a power-law stage to a de Sitter 
stage. The Chaplygin gas and its extensions were the 
unified models that have been more studied in the liter- 
ature [§j]-[l5j]. Nevertheless, preliminary results indicate 
that observationally the equation of state of dark energy 
still can not be determined precisely This information 
has produced several generalizations of those models, as 
for instance, the variable modified Chaplygin gas model 
with an equation of state, depending explicitly on the 
scale factor (ilj-pljl. 

At the present, there exist some controversy between 
interacting and unified models that we can express in the 
following question, are unified models more probable and 
satisfactory than interacting models profusely studied in 
the literature? Perhaps, this controversy does not exist 
and one could expect some sort of resemblance between 
these different models. In this case it would be partic- 
ularly interesting to find a relation between interacting 
and unified models. As far as we know, no study has been 



made in this direction. In the following we will examine, 
from the dynamical point of view, when interacting and 
unified models may be considered as similar ones. 

The paper is organized as follows: in section II we 
consider two fluids, dark matter and dark energy, with 
energy transfer, develop an effective one-fluid description 
and find the source equation for the total energy density. 
Then, we write the evolution equation for the effective 
barotropic index, introduce a separable interaction and 
investigate the conditions of stability for constant solu- 
tions. In section III we introduce the "linear interaction 
I" and describe some interacting models particularly sim- 
ple. After that we focus on the "linear interaction II" and 
find the exact scale factor and the effective barotropic 
index. Finally, we analyze a "general linear interaction" 
which induces a generalized ACDM model. In section 
IV we consider a "nonlinear interaction" which includes 
a rational function of the energy densities of both dark 
components homogeneous of degree one and separate the 
analysis into two main cases. There, we show that the ef- 
fective equation of state of two fluids with energy transfer 
includes the equations of state of the several generaliza- 
tions given for the Chaplygin gas. In section V we give a 
prescription to obtain an interacting model from a unified 
one. Finally, in section VI the conclusions are stated. 

II. DARK SECTOR EVOLUTION 

A. Effective one-fluid description 

Let us consider an expanding universe modeled by a 
mixture of two interacting fluids, namely dark matter and 
dark energy with energy densities p c and p x , and pres- 
sures p c and p x respectively. Due to the energy transfer 
between both dark components, they do not evolve sep- 
arately Q and the Einstein equations in a spatially flat 
FRW universe read 

3H 2 = p c + p Xl (1) 
Pc + Px + ZH(p c +p c + p x +p x ) = 0. (2) 

where a is the scale factor and H — a/ a. The conserva- 
tion equation ^ evidences the interaction among the 
components admitting the mutual exchange of energy 
and momentum. 

For the two dark components we assume equations of 
state p c = (7 C - l)p c and p x = (~/ x - l)p x , where the 
barotropic indices "f c and j x are constants. The dark 
matter is composed of nearly pressureless components 
with a barotropic index 7 C w 1 and the dark energy has a 
barotropic index satisfying the condition ^y x < j c . Many 
of our results will even be valid when we include the 
possibility of phantom dark energy 7 X < 0. The total 
energy density p and the conservation equation for the 
interacting two-fluid model are 

P = Pc+Px, (3) 
p' = -IcPc - JxPx, (4) 
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where the prime indicates differentiation with respect 
to the new time variable ' = djdr\ = d/3Hdt = 
d/d\n (a/ao) 3 and do is some value of reference for the 
scale factor. Solving the system of equations ©-(HI) we 
get the energy density of each dark component as a func- 
tion of p and its derivative p' 



7xp- 



A7 



Px 



IcP- 



A7 



(5) 



where A7 = 7 C — 7^ is the determinant of the linear 
equation system (J3j- ((4j) , being positive for our model. 
In turn, the ratio of the energies densities becomes r = 

Pc/Px = ~{lxP + P')/{lcP + (/)■ 

At this point, we introduce an energy transfer between 
the two fluids by separating the conservation equation 
for the system (0| into the two equations 



Pc 
I 

Px 



IcPc = -Q, 
- IxPx = Q- 



(6) 
(7) 



Here, we have consider a coupling with a factorized H 
dependence 3HQ where the interaction term Q, with 
dimensions of an energy density, generates the energy 
transfer between the two fluids. With this assumption 
the dynamics of p c and p x is dictated by the scale fac- 
tor instead of H . Differentiating the first or the second 
Eq. ([5]) and combining with the Eq. ((6]) or with the Eq. 
([7| , we obtain a second order differential equation for the 
total energy density 



p" + ilc + lx)p' + IdxP = QA7. 



(8) 



A similar equation was reported in [22j for the particular 
interaction Q = c\p c + C2p x , with c\ and c 2 constants. 

The interacting two-fluid model has been reduced to 
an effective one-fluid model with total energy density p 
and total pressure p = p c + p x , whose effective equation 
of state is 



P(P,P') 



-P- P 



(9) 



From the above point of view and the conservation of 
the total energy-momentum tensor of the system, we as- 
sume an effective one-fluid description with equation of 
state p = (7 — l)p, where the effective barotropic in- 
dex 7 = (7 c p c + r ) x Px)/p ranges between j x < 7 < j c . 
The effective conservation equation becomes p' + r yp = 0, 
hence the expression for the energy density of the dark 
components ([3]) are replaced by 

Pc = T—P' Px = — Pi ( 10 ) 



A7 



A7 



and the ratio becomes r = (7 — J x )/ (7c — 7)- Expressing 
these equations in terms of the energy density parameters 



7~7x 



A7 



7c - 7 
A 7 ' 



(11) 



In other words, given an interaction Q, the total energy 
density p of the effective one-fluid model is determined 
by solving the source equation Eq. ([8]). Once we know 
p, we are able to find the effective equation of state from 
Eq. ([9]) and the scale factor by integrating the Fried- 
mann equation 2>H 2 = p, without knowing p c and p x 
separately. Both energy densities are easily calculated 
by replacing p and p' into the Eq. ([5])- For instance, 
in the no interaction case, Q = 0, the energy density of 
the effective one-fluid model is p = Ci/a 37c + c 2 /a 37x , 
being p c the first term and p x the second one. Likewise 
p = (7c~ l)Pc + [lx - l)px = —p - p' is the effective 
equation of state. For any value of the constants c\ and 
C2 the total energy density p — > C2 /a 37:c , the scale factor 
a —> t 2 l Zlx and the power-law solution t 2 / 3 ^ becomes an 
attractor. Throughout the paper c;, c^, .... and bi, b[ .... 
with i = 0, 1, 2, 3...., will represent constants. 

Basically, we have shown that an interacting two-fluid 
model can be seen as an effective one-fluid model or 
equivalently considered as a unified one. Concerning this 
result, can the interacting two- fluid model or its unified 
version be derived from a Lagrangian? We observe that 
the dynamics of the unified or effective model is given by 
the two independent Einstein equations 



3H 2 = p, 



3H(p + p) = 0. 



(12) 



we get r = Q c /Q x . 



These equations cannot determine the three quantities a, 
p, and p because we have one degree of freedom. Usually, 
the system of equations (|12[) is closed with an equation 
of state p ~ p(p). When we assume that the effective 
energy-momentum tensor splits into two dark com- 
ponents, T ik = T9 k + Tg, the Eqs. ^ become Eqs. © 
and ©. The latter equations cannot determine the five 
quantities a, p Cl p x , p c and p x . To preserve the one degree 
of freedom of the unified model we have introduced an 
equation of state for each dark component p c = ("f c — l)p c 
and p x — (jx ~ l)p x - Then, by changing these equations 
of state we obtain a very large set of interacting models 
which are equivalent to a unified one, meaning that the 
decomposition into dark matter and dark energy is not 
unique. We should take into account this degeneration 
and after that, we should go to the central point of de- 
riving the interacting models from a Lagrangian. Some 
effort in this direction was reported in [23j . 

Although in unified cosmologies we have not dealt with 
the evolution of large-scale inhomogeneities, it is impor- 
tant to mention something about its non-homogeneous 
generalizations. This is an issue of interest because can- 
didates for the dark matter unification will only be valid 
if they ensure that initial perturbations can evolve into a 
deeply nonlinear regime to form a gravitational conden- 
sate of super-particles that can act like cold dark matter. 
In this sense one could follow the covariant and suffi- 
ciently general Zeldovich-like nonperturbative approach 
given in [8| . So that, it would be interesting to investigate 
in future works whether the equivalence between coupled 
and unified models holds also at the level of perturbation 
theory 
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B. Asymptotic stability 

The knowledge of stable solutions for the interacting 
two-fluid model is very useful because these attractor so- 
lutions determine the asymptotic behavior of the energy 
density of each dark component, so that, their contri- 
butions to the total energy density become constants. 
These scaling solutions are characterized by constant en- 
ergy density parameters fi cs , Q xs and they are reached 
for a broad range of initial conditions, hence alleviating 
the problem of coincidence. Other point of view consists 
in to look for a dynamical solution of the problem of coin- 
cidence such that the universe approaches to a stationary 
stage i.e., existence of attractor solutions for the evolu- 
tion equation of the effective barotropic index. Thus, on 
the attractor, the ratio r turns asymptotically constant, 

r s = Pcs/Pxs = Qcs/^xs = (7s - 7x)/(7c - 7s)- Consis- 
tently with both points of view the effective barotropic 
index tends to the asymptotic constant value j a , as we 
can see from 



7s 7c^C3 ~l~ 7a; ^xs 



Tslc + lx 

1 + r. 



(13) 



From these two relations we can extract two different 
types of solutions (i) for -f s ^ 0, we integrate the 
barotropic index 7 S = —2H/3H 2 and obtain the power- 
law expansion a = i 2 / 37s , (ii) for j s = 0, we have a 
final de Sitter stage, H = const, with £l cs =Q xs = and 
r s = --fxhe- 

As Eq. (jTB")) relates the constant density parameters 
i7 cs , f2 xs , the ratio r s and the effective barotropic in- 
dex 7 S , we will investigate the stability of the constant 
solution 7 S from the evolution equation of the effective 
barotropic index. Also, we will find the conditions of sta- 
bility for the solutions 7s and get the attractor when the 
energy transfer between both dark components is gener- 
ated by a separable interaction. To this end, we deduce 
the differential equation for 7 by differentiating p' = —7/) 
and by replacing p' and p" = ( r ) 2 — j')p into the source 
equation ([5]) so, we have 



7' - (7 - 7c)(7 - 7x) 



A7 
P 



(14) 



Firstly, we assume that a constant solution 7 = 7 S of 
the Eq. (fT4|) with j x < 7 S < 7e exists, and after that, we 
will impose the condition of stability so that 7 S be sta- 
ble. An interaction satisfying the existence requirement 
belongs to the class 



Q(7s) = (7 *-^ ( 7 " ) 



(15) 



with Q("f s ) < 0. The negative value of Q{j s ) indicates 
that the energy is being transferred from dark energy 
to dark matter, meaning that the latter component will 
dilute more slowly compared to its conserved evolution, 
p c oc a _37c , whereas the accelerated expansion of the uni- 
verse decreases compared with the noninteracting case, 



j x < j s . Interestingly enough, from Eq. (fj"3)) the result 
Q(~fs) < guarantees that the ratio r asymptotically 
tends to the constant value r s , thus alleviating the prob- 



lem of coincidence |24|-|28|. For the class of interactions 
(fT5j) , 7 S is a stationary solution of (fl"4)) and we obtain the 
scale factor a = i 2 / 37s by integrating 7 S = —2H/3H 2 . 

In what follows the analysis of stability will be 
restricted to interactions that have the form Q = 
Q(p c ,p x ,p' c ,p' x ,p,p',p"). By using the Eq. ^ and 
p' = —^ip with p" = (7 2 — 7')/0, we obtain that p CiX = 
Pc,x(p,p') and p' c x = p'c.x(p' t p")- then the interaction 
becomes Q = Q(7, 7', p). For simplicity we adopt sep- 
arability of Q, that is, Q = Q(7,7',p) = pQ(l,j') and 
write 



Q(i,i',p) 



(7 - 7c) (7 - lx) 



F(l,l')p, 



(16) 



where the function F depends on 7 and 7'. Several inter 
acting models analyzed in the literature are described b 
the interaction term ([To]) , see for instance [2 



351, \3 



For later proposals the separable Q also can be rewritten 
as 



Q(7,7» 



(17) 



by using the Eq. (|T0|) . In particular, for the F = const. 
case, the latter interaction ([TT]) will be investigated in 
detail in section IVB and there, we will show that the 
equation of state of the effective one-fluid model gener- 
ates several versions of the modified Chaplygin gas equa- 
tion of state. 

By combining Eqs. (|14j) and (|16j) we rewrite the evo- 
lution equation of the effective barotropic index as 



7 



"(7 -7c)(7- 7x)(F-l). 



(18) 



In conclusion, when the function F fulfills the two con- 
ditions 



and 



F(7 = 7„Y = 0) = 1, 



(7s ~ 7c) (7s ~ 7^)^(7^ 0) 
l + (7s-7c)(7s-7*)M7s,0) 



(19) 



< 0, 



(20) 

where F 1 and -Fy stand for the partial derivatives of F 
with respect to 7 and 7' respectively, then 7 S is a stable 
solution or an attractor. In other words, when the con- 
dition P^|) is satisfied the constant solution 7 S becomes 
a stationary solution of the Eq. (|18p. Besides, 7 S is sta- 
ble whenever the condition of stability ([2D")) is fulfilled. 
Hence the scaling solutions are attractors. 

When the interaction Q does not satisfies the existence 
requirement of belonging to the class (fl~5j) there is no con- 
stant solution of the evolution equation of 7. However, 
in section IV we will examine the structural stability for 
some interactions which do not satisfy the existence re- 
quirement, as for instance, the inhomogeneous nonlinear 
interaction. 
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III. LINEAR INTERACTION 

Various cosmological models investigated in the liter- 
ature are described by an interaction depending linearly 
on the energy densities p c , p x and the total energy density 
p [22j]- [iiil. From the beginning we have seen in Eq. ([5]), 
that p c and p x are linear functions of p and its derivative 
p'. This encourage us to investigate linear combinations 
of p c , p x , p and p' , namely the "linear interaction I" (Qi). 
We will study simple examples where the dark compo- 
nents interact with each other successively by only each 
term of Qi and review some of the models investigated 
with these particular couplings. For more generality, we 
will extend this study to the case of considering the "lin- 
ear interaction II" (Ql)- It has new terms proportionals 
to the first derivative of the energy densities p' c , p' x and 
a term proportional to the second derivative of the total 
energy density p" . It is motivated for the fact that the 
source equation (|8|) becomes a linear second order differ- 
ential equation for the interaction Ql- In this extended 
case we will solve the Friedmann equation and find the 
exact scalar factor together with the effective barotropic 
index. Also, we consider the possibility of adding a con- 
stant term to Ql and introduce the "general linear in- 
teraction" {Qql) with the intention of having a de Sitter 
scenario. After that we consider the contribution of a 
function of the scale factor and present the "more gen- 
eral linear interaction" {Q m gL) to obtain other different 
final stages. Below, we will investigate the outstanding 
aspects of those interactions. 



A. Linear interaction I 

Here we consider a linear combination of p c , p x , p and 
p' and define the "linear interaction I" in the following 
convenient form 

n , I \ (is ~7c) (7s -lx) . i > 
Ql(Pc,Px:P:P ) = Cl T P+C 2 ('J S ~ 7c)Pc 



When we impose the condition (jT9|) to the function 
, we obtain 



A7 



/ s (is - 7c) 7s - 7x) , , 01 , 

~ C 3 (ls ~ lx)Px - C 4 T P- (21) 

Is^l 

It has been obtained from Eq. (fT6"| by choosing the func- 
tion F, 

nw \ (ls-lc)(ls~lx) . Is -7c 
F;(7) = Cl— —yj- — r + C 2 - 



(l-lc)d-lx) 7 ~7c 



7s ~ lx , 7(7s ~ 7c) (7s - lx) , 00 ^ 

C 3 hC 4 -, r-, r-. 22) 

l~lx 7s(7-7c)(7-7x) 



Particular cases of the Qi (|21[) with ci = C4 = were 
previously inv estig ated in Refs. 22 1, [2!|-[30j], with c 2 = 
C3 = C4 = in 31|-[33[, with Ci = C3 = C4 = in [34| and 
with ci = c 2 = c 4 = in B-pcl. 



Cl + c 2 + c 3 + c 4 = 1. 



(23) 



Then, a constant solution 7 S of the Eq. (fT^|) exists pro- 
vided that the Eq. (|23|) be satisfied. Also, by combining 
Eqs. ©, ([21]) and ([23"]). we reduce the interaction ([3TJ to 
a function which depends only on the total energy density 
and its first derivative 

n , a hp + 7s' 1 i h - (is - 7c) (7s - 7*)] p' , 0/n 
Ql(P,P) = T ' ( 24 ) 



A7 

h = Ci(7 s - 7c) (7s - lx) ~ C2lx(ls - 7c) 
-C37c(7s ~lx)- 



(25) 



Inserting this Qi into the Eq. (HU, we find the two con- 
stant solutions 



li =7s, 



7j 



7c7x 



7s 



(26) 



In turn, the condition of stability (|20[) gives 

7s - 7, + < 0, (27) 

so, the constant solution 7 S is asymptotically stable pro- 
vided that 7,5 < 7 Z + or 7^ < idx — h. For positive energy 
densities also must be satisfied that j x < -f s < 7 Z + < -f c . 
Besides, for given values of 7 C and j x , the inequality (|2"7) 
bounds the range of constants c\ . c 2 , 03, C4 and the form 
of the Qi (|2"Tj) , yielding a stable cosmological model with 
the power-law expansion a = t 2 ^ 3js . 

By writing the source equation §E§ for the Qi (fM| . we 
have that 

+ 77 1 (7s 2 + Idx -h)p[ + ( lclx - h) Pl = 0, (28) 
and its general solution is 

pi = b ia - 3 ^ +b 2 a- 3 ^. (29) 

This total energy density has a vanishing limit for an 
expanding universe. Then the interacting model is finally 
realized when the general solution (|29[) is inserted into 
the energy density of each dark component ([5]) and the 
effective equation of state © 



Pel 



(7s- 


-7x)fcia" 


— + ( 7 + 


-yx)b 2 a 


-3,+ 


(30) 






A7 




* 


(7c- 


-7s)bia" 


- 3 ^+(7e- 


-7+)b 2 a 




(31) 






A7 




3 


(7s- 


- 1)« - 


f (it 


7s)&20" 




(32) 



The coupling between the two dark components mod- 
ifies dynamically, typical characteristics of p c and p x . 
In fact, the universe begins with a mix of dark mat- 
ter (|30p and dark energy (|31[) represented approximately 
by the unstable energy densities p c oc (7^ — 7a:)a _37!+ 
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and p x oc (jy c — 7 ( + )a 37 ; respectively. After that, the 
instability of the constant solution 7 ; + induces the uni- 
verse to evolve from that unstable era, characterized by 
r + = (7+ — 7i) /(7c — 7 ; + ), to a stable final stage where 
the dark matter and dark energy densities are domi- 
nated by the stable components p c oc (7 S — 7:r )a~ 37s and 
Px °c (7c — 7s) a ~ 37s - The stable solution 7 S is associated 
to an asymptotically stable ratio r s = (jy a — i x )/ (lc — 7s) 
with the expected result r + > r s , showing that the lin- 
ear interaction alleviates the problem of coincidence. In 
turn, the scale factor interpolates between the unstable 
stage, evolving as a oc i 2 / 37 i , and the stable stage evolv- 
ing as a oc i 2 / 37s . Meanwhile the effective equation of 
state (|32[) plays the role of a peculiar fluid at the initial 
stage, while at late times it turns into the equation of 
state of a perfect fluid. For large scale factors the quan- 
tities p, p c , p x , p' , p and Qi behave as a _37s , in this way, 
the original evolution equations ([6]) and (|7|) for the dark 
components become algebraic equations. 



B. Linear examples 

Now, we analyze four simple cases by considering sep- 
arately, each term of the Qi (|2ip . In all these examples 
we select the four function ^(7) in such a way that the 
condition (fT9")l is satisfied identically 



1. c 2 = c 3 = c 4 = 



Q, 



(7s ~ 7c) (is ~ lx) 

A7 



{ls-lc){ls-lx) 

(7-7c)(7-7a) 



(33) 

This negative interaction Q p < represents an energy 
transfer from dark energy to dark matter. The inter- 
action Q p between quintessence and a pressureless com- 
ponent in a spatially flat FRW cosmology produces a 
transition from a phase dominated by dark matter to an 
accelerated expansion phase dominated by dark energy 
[3l | -[33 j . Simultaneously the interaction Q p alleviates 
the problem of coincidence of the present universe Q . 
By imposing the condition (|20p on the function F p 



1 



1 



7s - 7c 



<0, 



7s - lx 



(34) 



we find that the solution 7 = 7 S is an attractor provided 
lx < Is < (7c + 7x)/2 < 7c- Inserting Q p into the source 
equation (jHJ we obtain the total energy density and the 
effective equation of state (0) 

p = M~ 37s + & 2 a" 3(7c+7j: " 7s) , (35) 
p = ( 7s - l)p + ( 7c + lx - 2 7s )& 2 a- 3(7c+7 *- 7s) • (36) 

For any value of the initial conditions &i,&2 and large 
scale factor, the total energy density p — > bi/a 3ls and 
the effective equation of state (|3"6"|) adopts the barotropic 
perfect fluid form p tn ( 7s — l)p with a — s- t 2 / 37s . Hence, 
the power-law expansion is asymptotically stable. 



The dark matter and dark energy densities ([5]) can be 
written as 

pc = (7s-7x)fria- 3 ^+(7^-7 s )fr2n- 3 <^+-^-T s ) ^ 

Thus, the ratio r p tends to r s = {l s — l x )/{lc~ls)i being 
r s an attractor. 



2. ci 



C3 



c 4 = 

= (7s - 7c)Pc 



7s -7c 



(39) 

7 -7c 

For this interaction, with Q Pc < [34| , the condition (|2T))) 
is not satisfied 



<(7s) = (7c-7s)" 1 >0, 



(40) 



and the power-law solution a = £ 2 / 37s , with 7^ restricted 
to the interval ^ x < j s < 7 C , is unstable. This model 
contains serious instabilities on the perturbations of the 
dark energy component [35| . 

In this example the solution of the source equation ([8j) 
with the interaction Q Pc and the effective equation of 
state El) are 



p = b ia - 3 ^ + b 2 a- 3 ^ , 

P = (7s - l)p + (7* - ls)b2a~ 31 * 



(41) 
(42) 



For any value of the initial conditions b\ ,62 and for large 
scale factor, the total energy density (|4T|) has the limit 
p Pc — > 62/a 37x , meaning that a — » f 2 / 3 ^ since j x < 7s . 
Accordingly, the effective equation of state (|42j) becomes 
that of the dark energy p sa (7^ — l)p indicating that the 
interaction Q Pc would not be adequate to describe the 
evolution of dark components. The model seems to be 
completely dominated by the dark energy. In fact, the 
energy densities E)) are 



Pc 
Px 



(is -7z)M 373 
A7 

(7c-7s)M~ 37s 



A7 



+ b 2 a~ 3 ^ 



(43) 
(44) 



and the ratio r„ 



oc a 



-3(75-7*) 



0. Then, at late times, 



the interacting two-fluid model with energy transfer Q Pa 
would not solve the problem of coincidence and it would 
not be suitable for to fit the present observations. How- 
ever, this coupling can work when it is combined linearly 
with some of the other parts of Qi . 



3. c\ = c 2 = C4 = 

Qp, = -(is ~ lx)Px 



7s -lx 

l~lx'' 



(45) 



with Q p < 0. This interaction was examined in several 
papers 24 1- [2^]. Now the condition of stability ([2"U]l is 
satisfied 



^(LCTs) = (lx — 7s)" 1 < 0, 



(46) 
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and the solution 7 S is stable. By solving the source equa- 
tion (|8) for Q Px and using the Eq. we obtain the 
total energy density and the effective equation of state, 
they read 



P 
P 



(is ~ l)p+ (7c - 7s 



j 2 a 



-37 c 



(47) 
(48) 



For any value of the initial conditions bi,b 2 and large scale 
factor, the total energy density p — > ci/a 37s and a — > 
t 2 /37 3 because 7 S < 7 C . The effective equation of state 
behaves as p w (7^ — l)p showing that the interacting 
model is dominated by the attractor 7 S . 

The dark matter and dark energy densities ([5]) are 
given by 



(is -lx)bia 373 _ 3 
Pc = — 1- b 2 a J7c , 



Px 



A7 

(7c - Is) M~ 37g 

A7 



(49) 
(50) 



showing that the ratio r Px — > r s = (i s — l x )/(lc — 7s) 
on the attractor. Then, the interaction Q P:c may rep- 
resent adequately a coupled model of dark matter and 
dark energy. In fact, the ratio of these components has 
enough parameters to be adapted to the observations, 
consequently this interacting model may be a candidate 
to alleviate the problem of coincidence. A cosmolog- 
ical model with the above characteristic was proposed 
for the current universe which consists of noninteracting 
baryonic matter and interacting dark components [24| . 
There it was used two interacting fluids in the dark sec- 
tor with constant barotropic indices. The energy trans- 
fer was taken proportional to the dark energy density 
Q Px and it was shown that the model leads to a correct 
behavior which is expected for a viable scenario of the 
present universe, as for instance, the deceleration param- 
eter, density parameters and luminosity distance. Also 
the interaction p Px was used to show that the overall en- 
ergy transfer should go from dark energy to dark matter 
if the second law of thermodynamics and Le Chatelier- 
Braun principle are to be fulfilled, guaranteeing that the 
ratio r Px asymptotically tends to a constant, thus allevi- 
ating the problem of coincidence [26| . The evolution of a 
viscous cosmology model was also analyzed by employing 
an energy transfer between the dark components induced 
by the interaction Q Px [27| . 



4. ci = c 2 = c 3 = 

q / = (lc - 7s)(7s - lx) , p i _ l(ls-lc)(ls-lx) 

9 7sA7 ' p 7s(7-7c)(7-7*) ' 

(51) 

By using p' = —jp < 0, we see that the interaction Q p i is 
negative. When the condition of stability (|2"0")) is imposed 
to the function F p i , we get 



7c7x - 7; 



7s(7s -7c) (7s -lx) 



<0, 



(52) 



which means that 7,^72; > 7^. This inequality determines 
the range of values of the attractor 7 S . Solving the source 
equation ((8]) for Q p i and inserting it into the Eq. ((9]), we 
find the total energy density and the effective equation 
of state 

p = M~ 37s + b 2 a-^^ h ' , (53) 
p = ( 7 s - l)p + 77 1 (7c7* ~ 7s 2 )M- 3W7 =. (54) 

When 7 C j x > 7^ is satisfied whatever be the initial 
conditions b±, b 2 the total energy density has the limit 
p — > bi/a? ls for large scale factor, evidencing that j s is 
an attractor. 

The dark matter and dark energy densities ([5]) are 
given by 

Pc — A7 ' \ 00 ) 

= (7c -7, )bia- 3 ^ +7c77 1 hs ~7x )b 2 g-^--<- ^ 



and r r 



(is — lx)/(lc — 7s) orL the attractor. As far 



as we know, this interacting model was not investigated 
in the literature. It appears as a feasible candidate to 
be considered for to describe the evolution of the dark 
components and to alleviate the problem of coincidence. 



C. Linear interaction II 

To enlarge the set of linear interactions, we take into 
account a coupling with the first derivative of dark mat- 
ter and dark energy densities, p' c and p' x . Coming back 
to Eq. we observe that the above assumption in- 

troduces a dependence with the second derivative of the 
total energy density. Then, we generalize the Qi adding 
it these new terms and define the " linear interaction II" , 
Ql, in such a way that it verifies the conditions of sta- 
bility (ri"9"]) - (|2"Ul) . As far as we know, this kind of coupling 
was not investigated in the literature. So that, we will 
analyze in detail the Ql. 

We start building the Ql by combining linearly the 
quantities p c , p x , p' c , p' x , p, p' and p" 

Ql = cip c + c 2 p x + c 3 p' c + C4p x + c 5 p + c 6 p' + c 7 p". (57) 

By using the Eqs. ([5]) the first and second terms of the Q l 
reduce to a linear combination of the total energy density 
and its derivative whereas the third and fourth terms 
reduce to a linear combination of the first and second 
derivatives of the total energy density. Then, rearranging 
all the terms in the Eq. (|57| . the Ql can be reduced to 
a linear combination of the basis elements p, p' and p". 
Finally we get 



c' 2 p' 



c' sP ", 



(58) 



where the constants c[ are linear combinations of the con 
stants Cj in the Ql ([BT|) . On the other hand, the Ql 
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also can be obtained from Eq. (fl6]) by selecting the func- 
tion F L (-f,Y) as 



Ft 



(is - 7c) (73 - 7s) 
(7-7c)(7-7r) 



u 1 u 7 . , 7 2 - i 
bi + o 2 h 03 ^ — 

Is It 



(59) 

where we have used that p' = — 7/5, p" = (7 2 — 7')^. The 
coefficients 



61 = 
b 2 = 
b 3 = 

satisfy the constrain 



A 7 ci 



(7s -7c) (7s -7*)' 
7s A7 4 
(7s -7c) (7« -7*)' 



7 S 2 A 7 4 



(7s -7c)(7s-7x)' 



(60) 
(61) 
(62) 

(63) 



b 1 + b 2 + b 3 = 1, 



after imposing the condition (|T9l) to the function (|59|) . 
By combining Eqs. ([58|) . ([59f and (|63|) . wc obtain the 
final Ql 

_up + 77 1 [w + u - (7 S - 7 C )(7 S - 7 X )] p' + v% 2 p" 



with 



A7 



u = (is - 7c) (7s - 1x)bi, 

V = {is - 7c)(7s ~ 7s)&3- 



(64) 

(65) 
(66) 



Coming back to the Eq. (fT4]) , the Ql generates the two 
constant solutions 



+ Idx ~ U 

1l =7s, 1 l =1s—t— 



(67) 



7| - v 

while the condition of stability (f2T)]) yields 

7s - 7 £ < 0, (68) 

with the additional requirement j x < is < it < 1c- 

For Ql the exact general solution of the source 

equation ([5]) is given by 



p L = b[a- 3 ^ + b' 2 a- 3 < . 



(69) 



Then, the Ql induces a stable cosmological model with 
a final behavior described by the power-law expansion, 
a = t 2 / 373 . In addition, the energy density of each dark 
component (JSJ) and the effective equation of state (|9]) are 

„ _ (7 3 ~7,)fc' 1 q- 3 ^+(7+-7,)^q' 3 ^ (m) 

PcL — A7 , I. ' u; 

_ ( 7e - 7a ) ' 1 a"^+( 7e - 7 +)b^- 3 "'L 
Ae£ — A 7 ' V / 

PL = (is - I) PL + (it - 1s)b' 2 a^l . (72) 



For large scale factors, similarly to the Qi case, the quan- 
tities p, p c , p x , p' c , p' x , p', p", p and the Ql behave as 
a _37s , in this way, the original evolution equations of the 
dark components and ([7]) become algebraic equations. 
In the initial regimen the above quantities behave as in 
the Qi case with 7 ; + substituted by 7^. 

Finally, the Friedmann equation 3H 2 = pl for the 
source (|69l) is implicitly solved and one finds the scale 
factor 



a-L 
t = 



[u> sinhAr 
2 



2/3(-yt-7a) 



(73) 



^2 (it 



Is) 



J [a;sinhAr] 7s/(7 ^" 7s) dr,(74) 



where uo 2 = b^/b^, see details in [3(|. Due to 7^ > is, 
the latter equation shows that the variables t and r have 
the same asymptotic limits. Then, it is appropriate to 
investigate the scale factor and the remaining quantities 
in the two asymptotic regimes. The effective barotropic 
index reads 



1L 



it 



- 7 S sinh luAt 



cosh wAt 



(75) 



so as t grows the model interpolates between the initial 
7^ and the final 7 S values. Eqs. ([75|) -([75 |) allows us to ex- 
press the dark matter and dark energy densities (fTDj) . the 
total energy density (|69|) . the ratio tl and the effective 
pressure © as functions of the new time r. In particu- 
lar, at early and later times, the asymptotic limits of the 
ratio tl become 



1l -7s 



7c 



it 



7s ~7s 
7c - 7s ' 



(76) 



These ratios satisfy the crucial relation > r s , thus 
the Ql gives the possibility of alleviating the problem of 
coincidence. 



D. General linear interaction and ACDM model 

We complete the subject of linear interaction by en- 
larging the basis elements with a constant, so that, the 
new base will be c, p, p', p". After that we generalize 
this basis elements by introducing a function of the scale 
factor instead of the constant c. Although, the effective 
one-fluid model is able to mimic the essential features of 
the ACDM cosmological model, clearly, the introduction 
of both modifications could produce radical alternatives 
to the ACDM model. 

With this aim in mind, we first introduce the " general 
linear interaction" , Q g L , by adding the constant Qq / A7 
to the Ql, so 



QgL - 



Qi 



(77) 



where we have assumed that the constrain (|63p holds for 
the constants in the Ql- Obviously the condition of sta- 
bility docs not hold. A particular type of the linear 
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combination (|77|). Q = cq + c\p c + c 2 p x , was analyzed in 
[37| . Here, we investigate the consequences of the Q 9 l on 
the ACDM model. Combining Eqs. flEJ, ((53)1. (f67]> and 
((77)) the source equation ((5) for the total energy density 
becomes 



(1 - «7 S 2 )p" + 7 S 1 (is + Idx 



v p 



+ ildx - u)p = Q . 



(78) 



Its general solution and the corresponding effective equa- 
tion of state (©, are 



PgL = A e// + b' ia -^ + b' 2 a-^t, 



(79) 



PgL 



-A e//7s + ( 7s - l) PgL + ( 7 + - 7s )^a" 37i , (80) 



where A e // = Qols/ltils ~ v ) ^ s the effective cosmologi- 
cal constant, induced by the constant term Qo/ A 7 in the 
Q ff L- In order to A e // > 0, we need to choose the param- 
eters of the interacting model such that Qo/(l 2 —v)>0. 

At late times the total energy density has the limit 
p g L — > ^-eff and the effective equation of state becomes 
p w —A e ff. Thus, the effective one-fluid model can be 
associated with a unified dark sector model whose scale 
factor interpolates between a power-law phase and a de 
Sitter stage H = Ho = y/3/A e ff. Being Ho the attractor 
solution in the phase space for any value of the initial 
conditions b[ and b' 2 - 

From Eq. ([5]) , the dark matter and dark energy densi- 
ties p C L and p X L are 



PcL 



7*A 



eff 



A7 



PcL, 



PxL 



7cA e// 
A7 



PxL- (81) 



It is interesting to comment that, in general, when the 
total energy density tends asymptotically to A e ff for 
a — >• 00, the Eq. ([5]) shows that the dark matter en- 
ergy density has the final limit p c — > —j x A e ff/A'y < 0. 
To relieve this problem we may assume a phantom equa- 
tion state for the dark energy, with "f x < 0, so the energy 
densities ((8T)) become positive and the ratio r g L tends to 
= — 7x/ 7 c > 0. Curiously, the Q s l dresses the bare 
phantom dark energy density p x and gives a decreasing 
dark energy density ((81)) with the final limit 7 c A e ///A 7 . 
In the cases which we admit that 7 S or 7T, or both, are 
negatives one or the two dark components have a final 
phantom phase with increasing energy densities. 

A "more general linear interaction", Q m gL, can be in- 
troduced by adding a well-behaved function f(j]) to the 
Q gL , so we get 



IragL 



lt(l 2 s- V ) 
ls&l 



A 



eff 



m 
A7 



Qi 



(82) 



The general solution of the source equation 
tained by finding the particular solution of (1— v~/ s 



is ob- 

; 2 )p"+ 

Is 1 {is + Idx -u-v) p' + (7 c7k - u)p = f(r)) and 
adding it to the total energy density (|79[) . Essentially, 



for large scale factors, the behavior of the effective one- 
fluid model is defined by the relative weight between the 
constant 1l{1s ~ v )^ef / / 7s A 7 and the function f(ij). 
In this way, the additional / term determines the final 
behavior of the more general interaction model and how 
much it deviates from the cosmological ACDM model. 
The case where the coupling (|82p reduces to the / term 
was investigated in [38| . 



IV. NONLINEAR INTERACTION 

Let us assume that the energy transfer between the 
dark matter and dark energy components is produced by 
the following "nonlinear interaction", 



SnL 



c 8 p c + c 9 p c p x 



ClOPx 



Ql 



f(v)p" 

A 7 : 



(83) 



where the Ql is giving by Eq. (|64p. /(^)/0 !y is a nonlin- 
ear atypical term proportional to a well-behaved function 
/(??), which depends on the scale factor 77 = In a 3 , and 
to p" '. The constant v will be determined later on so 
that the source equation (jSJ) can be recast in a solvable 
equation. 

So far we have examined interactions depending lin- 
early on: the dark matter and dark energy densities, 
their first derivatives, the total energy density and its 
derivatives up to second order. However, it will be in- 
teresting to consider an interaction term which includes 
a rational function of the energy densities p c and p x ho- 
mogeneous of degree one, as for instance the first three 
terms in the RHS of the Eq. (|83|), It could be consid- 
ered as a possible generalization of the Qi ((!M)) and the 
Ql (|53)l . In subsection IVB we will see that the rational 
interaction terms are of particular importance. Besides, 
a coupling depending on the scale factor will be useful so 
that the equations of state of the effective one- fluid model 
produces equations of state which extend and generalize 
those of the Chaplygin gas &- 11511 and the variable mod- 
ified Chaplygin gas model |16|-)21|. 

By using Eq. ((SJ) we have shown that the Ql reduces 
to a linear combination of the three elements of the base 
p = c^p + c'zp' + c'zp", (sec Eqs. ((SH)) and ([53])). Following 
similar arguments, the three terms in the numerator of 
the nonlinear part of the Eq. ([55]). c§p^ + Cgp c p x + c\op x , 
become a linear combination of p 2 , pp' and p' 2 . Rear- 
ranging all these terms, we obtain the final form of the 
nonlinear interaction 



n C{p' 2 + p (4p + dip? + c'jp") + fRA . 
QnL = p~A7f ' (84) 

where the c" are linear combinations of the remaining 
constants. Below we will see that the nonlinear terms 
arising from the homogeneous function of degree one in 
the energy densities p c , p x and f(rj)p 1 ' produce original 
results, as for instance, the effective equation of state 
describes universes which have power-law and the Sitter 
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expansions. Also these terms are really the responsible 
that the effective equation of state has a reminiscence of 
the Chaplygin gas, including its different versions. From 
Eqs. © and (|84p . we obtain the nonlinear differential 
equation for p 



n , 7c + 7x-C<j , 
PP + 1 -J, PP 



Iclx - Cj 2 

1-c" 9 



1 



JJP' 



(85) 



In the particular case where the first and last terms in 
the RHS of the Q n L (|84j) vanish simultaneously, the Eq. 
(|85[) becomes a homogeneous linear differential equation 
for p. However, in other cases, the general solution of the 
source equation (|55|) will be obtained from a nonlinear 
superposition of the two basis solutions of a second order 
linear differential equation. 

Renaming the four constants in the Eq. (|85[) by 61, b 2 , 
63 and 64 respectively (from left to right) and changing 
to the new variable x = p( 1+f>2 ' with 62 7^ —1 or c" + 
c'l 7^ 1, the Eq. (|8"5j) turns into the equation of a forced 
dissipative (61 > 0) or antidissipativc (61 < 0) linear 
oscillator 



x" + bW + 63(1 + h)x = 64(1 + b 2 )f(v), 



(86) 



where we have chosen v = —b 2 = c![j{\ — c'l). If X\h 
and X2h are the two basis solutions of the homogeneous 
Eq. (|86|) . then the general solutions of the Eq. (|85|) can 
be written as a nonlinear superposition of these basis 
solutions p n L = (b^xih + bex 2 h + x p ) 1 ^ 1+b2 - ) . There x p 
is the particular solution of the Eq. (|56")l and 65, b$ are 
integration constants. From now on, the Eq. (|86|) will 
substitute the source equation ([5]) for the Q n L- 

In what follows we will study the solution of the Eq. 
(|86|) considering separately: the homogeneous case / = 0, 
the example of Chaplygin gas and the inhomogeneous 
case / 7^ 0. The latter it will be divided into two parts, 
the / = /o = const, case, and the / = f{a) case where 
the "relaxed Chaplygin gas model" will emerge. 



A. The homogeneous case 

In the / = case, the source equation Eq. (|86|) be- 
comes homogeneous and the Q n L reduces to 



Qh = 



c' lP - l p 12 + c' 2 p + c' 3 p' + c' iP " 
A7 



(87) 



By inserting the "homogeneous nonlinear interaction" 
(Qh) into the Eq. we easily get the general so- 

lution Xh — 65a 3A «^ + bQa 3X "L. Then, coming back to 
original variable ph = x h ^ 1+b2 \ we have the total energy 
density, it reads 



p h = b 5 a 3X ^L + b 6 a 



1/(1+62) 



(88) 



where b 5 and b 6 are integration constants while X nL and 
X^ L are the characteristic roots of the linear source equa- 
tion (ESI) 



Kil - 



(89) 



On the other hand, the dark matter and dark energy 
densities ([3]) are 



Pch = -D 
Pxh = D 



3A+ 

[A^ + 7*(l + fc)]p + &eAA 



[K L +lc{l + b 2 )]p + b^\ 

P 



(90) 
(91) 



where D =[(l + b 2 )A^ 1 and A A = X+ L - X~ L . 

We have started from an interacting two-fluid model 
and finally obtained an unified cosmological model where 
the dark matter and the dark energy evolve as an effective 
one-fluid. Now, combining Eqs. © and we find the 
effective equation of state 



Ph 



1 



a; 



1 + 62 



Ph 



b^AX a 



3A" 1 



(1 



(92) 



Depending on the values of the parameters b\, b 2 , 
and &3 the total energy density (|55|) behaves asymptot- 
ically as ph — > a 3A "t^ 1+b2 - ) in the limit of large scale 
factors, meaning that a -t t~ 2( - 1+b2 ^ 3X r,L _ p or fr 5 = 
(+) or & 6 = (-), the model includes the exact power- 
law expansions a ± = i~ 2 ( 1+b2 )/ 3A n.L . In particular, when 
(1 + b 2 )/3X^ L > 0, we have a final phantom phase en- 
tailing that the scale factor expands so quickly that the 
scalar curvature R — > 00 in the limit a — > 00 and it 
reaches a = 00 in a finite amount of proper time. In the 
special cases that X^ L = and A^ L < 0, namely b\ > 
and 63 = or bi < and 63 = (sec Eq. (JHSJ)), we 
have a final de Sitter stage with an effective cosmolog- 
ical constant given by the limit ph — > A e // = &g + 



Or p h -> A e ff 



.1/(1+62) 



These models include the 



modified Chaplygin gas introduced in Ref. [Ill ] where it 
was proposed the equation of state p = Ap — B/p n , with 
n > 1, and the parameters A and B were constrained 
to be positive. When both, X^ L ^ and X~ L ^ 0, the 
equation of state (|92|) contains those which characterize 
various of the variable modified Chaplygin gas models 
investigated in Refs. [l6j]-[2l[. in the next subsection 
we concentrate on a specific nonlinear interaction and 
study the two associated cases where the characteristic 
root A+ L vanishes. 



B. Nonlinear examples 

Here we investigate the interesting case where Qh takes 
the form 



, PcPx 

P 



(93) 
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with a constant, see [6J and references therein. The in- 
teraction (|93p also can be obtained from the Eq. (fTTj) 
by choosing the constant function F(jy) = —a^c/A"/. 
From Eqs. ([55)1 and ([9"3")) we identify the three coefficients 
h = (7c+7x)(l + ^2), b 2 = aj c /Aj and 63 = 707^(1 + 62). 
Then, the characteristic roots ([89]) of the source equation 
(Ea read 



A^ = -7*(l + fe). 



A 



-7c(l+&2). (94) 



Introducing these roots into the homogeneous general so- 
lution ([88)) and the effective equation of state (1921) . they 
reduce to 



6 -3 7o (l+6 2 ) +6 -37.(1+6,) 



1/(1+62) 



P = (7c - 1)P - &6A7- 



-3 7ic (l+6 2 ) 
p6 2 



(95) 
(96) 



This effective equation of state can be identified with 
those which were used to build variable modified Chap- 
lygin gas models [l6j-f2l|. 

1. modified Chaplygin gas 

We make an adequate selection of the parameters so 
that the interaction ([9"3")) be focused on a dark energy 
component described by some kind of vacuum energy 
density i.e., j x = 0. Then, b 3 = 7 c 7 x (l+&2)=0, A7 = j c , 
62 = ol 7^ — 1) ^nL = 0, and the second term in Eq. 
([53)) becomes constant. Assuming that the dark mat- 
ter component is nearly pressureless, we may associate 
it to a barotropic fluid with a free constant parameter 
7 C = 70 ~ 1. Then the energy density (|95p abbreviates 
to 



B /a \ 3 To(i+") 



7o 



\ a / 



1/1+c 



(97) 



where the new constants B and eto are redefinitions of the 
old integration constants 65 and be. Hence, by replacing 
this energy density in Eq. we obtain the equation of 
state of the one effective fluid 



P = (7o - 1)P 



B 



(98) 



It characterizes several unified cosmologies implemented 
with Chaplygin gases as t he g eneralized, extended, modi- 
fied and enlarged ones [l]-[l5(. Also, these unified models 
along with some others ones [39[-|45|, were extensively 
used to describe unified versions of dark matter and dark 
energy. 

2. reduced unified model 

For 70 = 0, the effective one-fluid model described by 
the expressions (jpTJ) and (|98[) is not valid. In this "re- 
duced unified model" the source equation (|86|) must be 



solved again by making the particular choices b\ = 63 = 
and 62 = a. After comparing the latter particular 
choices with the coefficients of the nonlinear version of 
the source equation ([55)) we obtain the coefficients of 
the corresponding "reduced nonlinear interaction" (Q r ), 
C3 = 7 C + j x , c'2 = idx; / = and c'{ = a(c'l - 1) re- 
spectively Then, from Eqs. ([55)) and (|84|) . we get the 
final Q r 



Qr — 



a(c'l - \)p- x p 12 + 7c7xP + (7c + lx)p' + 



4p" 



A7 



(99) 

Now, the source equation ([go]) reduces to x" = and its 
general solution is x = b% + 6977. Thus, one finds that 
the total energy density p = x 1 ^ 1 " 1 "") has a logarithmic 
dependence with the scale factor 

l/(l+a) 



P = Po 



±l + 61n 



(100) 



where the constants po and b are related to the old inte- 
gration constants 6s and 69. In turn by using Eqs. ([5]) 
and @, we find the dark matter and dark energy den- 
sities, and the equation of state of the effective one-fluid 
model are 



Pc 



Ar 



P 



P 
A7 



1+Q 



P 
A7 

-P- 



7c + 



1 + Q' 

bpo ( Po 



l+a 



1 + a 



(101) 
(102) 
(103) 



The last equation of state also can be seen as a gener- 
alization of the polytropic equation of state p = Kp 7p , 
where if is a constant and 7 P is the polytropic index. On 
the other hand, by integrating the Friedmann equation 
with the total energy density (|100[) . we determine the 
exact scale factor 



a = a exp — 
3b 



±1 



2(1 + a) 



(104) 

where we have set t = 0. Note that the final results 
(|100p - (|104[) arc independent of c'{, then the interaction 
term p" dos not contribute to the evolution of this re- 
duced unified model. 



C. The inhomogeneous case and the "relaxed 
Chaplygin gas model" 

Here we are going to consider the " inhomogeneous non- 
linear interaction" (Qi), 



Qi 



c'lp-'p 12 + 4p + dip' + c'lp" + fiv)^ 1 -^ 



A7 



(105) 
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which has been obtained from the Eq. ([84)1 by choosing 
v = —b 2 = c"/ (1— c'l). The source equation (jg6")l becomes 
inhomogeneous and the general solution is obtained by 
adding a particular solution x p to the homogeneous so- 
lution Xh = b^a^riL + 6 6 a 3A "i , then the general solution 



Pi = (x h + xp) 1 /^^ of the Eq. 



+ Xr 



reads 
1/(1+62) 



(106) 



and the equation of state © takes the form 



Pi 



1 



A 



nL 



1 



6 fi AA 



,3A" 



a; 



(1 + W 

(107) 

(i) For f = fo = const, the particular solution x p — 
64/0/63 is constant and the total energy density (j 106[) 
along with the effective equation of state (|107[) describe 
a " double unified model" , in a sense that initially the uni- 
verse is dominated by the " two" terms inside the bracket 
of Pl ps (b 5 a 3X nL +6 6 a 3A ni) 1 /(i+62). These "two" terms 
can be seen as a " nonlinear mixture of two fluids" . But 
at late times the universe is dominated by a vacuum en- 
ergy pi ps Xp^ 1+ and has a de Sitter expansion. So, 
this (i) case includes various generalizations of the mod- 
ified Chaplygin gas model investigated in Refs. [l6|-(2l|. 
When the constants 65 or b$ vanishes the effective equa- 
tion of state (|107p turns in Eq. ([55)1 and the double uni- 
fied model produces different versions of the Chaplygin 

gas p-pii. 

(ii) For f(j]) ^ const, and r\ = In a, the particular 
solution x p and consequently the numerator of the last 
term in the RHS of the effective equation of state (|107p 
become arbitrary functions of the scale factor. Then, a 
fluid obeying the Eq. (| 107[) defines a " relaxed Chaplygin 
gas model" . 

Finally, we investigate the structural stability of the 
solutions corresponding to the source equation (|86p for a 
well-behaved function f(n), which will be considered as 
a degree of freedom in the problem of stability. 

We use the analogy with the classical potential prob- 
lem by writing the Eq. (|86[) as an equation of motion for 
a dissipativc mechanical system, namely, 



d 
dij 



where 



V{x,n) 



— + V(x, V ) 



63(1 + 62) 



= -D{x,x\n), 



_i f 2 

bV 



D(x,x',n) = b i {l + b 2 )f'x + b 1 x'' 



(108) 

(109) 
(110) 



V(x,T]) is the "potential" and D(x,x',n) is the "dissipa- 
tion" of the " equivalent mechanical system" . In addition, 
we assume that the potential V is inferiorly bounded by 
a function Vj,, it has a finite limit for rj — > 00 and the 
dissipation D is positive definite in the same limit. 



When 63(1 + 62) > 0, the potential becomes bounded 
by the function VJ,, thus V > where 



V h 



2b: 



(in) 



This bound is obtained by evaluating the potential (|109j) 
on the limit solution xu m = b^f /b% of the source equation 
||SBJ> for 77 -> 00, or p Um = (64//6 3 ) 1/(1+b2) . Therefore, 
for a dissipativc mechanical system with b\ > 0, the set 
of functions f{rf) such that / — > f for r\ — > 00 and 64(1 + 
62)/' > define a positive definite dissipation, D > 0. 
In this context, a solution of the source equation which 
have the limit p — > p m in may be considered as a " stable 
solution" . 



V. STARTING WITH UNIFIED MODELS AND 
ENDING WITH INTERACTING ONES 

Eqs. ©-(Ill) introduce an alternative interpretation of 
unified models by associating them with interacting mod- 
els. Several kinds of unified models were extensively in- 
vestigated in the literature [li|-[45j], for instance: the 
Chaplygin gas and purely kinetic, quintessence, k-essence 
and DBI cosmologies. For purely kinetic cosmologies we 
mean that the pressure and energy density of each field 
involve only the time derivative of the same field. 

At this point we consider two options, in the first one 
the equation of state of the unified model has the form 
p = p(p), as it is for the Chaplygin gas (see the Eq. 
(1981) ) or we have purely kinetic quintessence, k-cssence 
and DBI fields. These are described by the equations of 
state p q = p q - 2^, p k = -F{4> 2 {p k )) where T = F(4> 2 ) 
is the kinetic function and <fr is the k-essence field, and 
Pdbi = -Vq + (p- V a )/{pDBifQ_- V fo + 1) for Dirac- 
Born-Infcld (DBI) cosmologies [45(. The parameter Vq 
represents a constant potential for the three fields and fo 
is the constant warp factor of the metric. In all the above 
unified cosmological models, after integrating the conser- 
vation equation of each component, the total energy den- 
sity depends only on the scale factor, p = p(a). In the 
second option we include the cases where the equation 
of state p — p(p) is given but we do not know p = p(a) 
explicitly. 

To clarify the alternative interpretation, let us assume 
that we have found the energy density p = p(a) of the 
fluid which characterizes a unified model. Then the Eq. 
§5§ associates this unified model with a two-fluid model 
whose components have energy densities p c (a) and p x (a) 
respectively. Besides, inserting p(a) into the source equa- 
tion ©, we can calculate the interaction Q = Q(a). 
Therefore, a one-fluid model with energy density p(a) 
interpolating between dark matter and dark energy sce- 
narios, can be expressed as a two-fluid model with energy 
transfer. In other words they share the same Fricdmann 
equation and the same scale factor, so they define the 
same geometry. 



13 



In the case we know the equation of state p = p(p) 
of the unified model, instead of p(a), we substitute the 
former into Eq. ((5]) and obtain the energy density of each 
fluid 



Pc 



P+ (1 ~ 1x)P 

A 7 



Px = 



-P + (7c - 1)P 

A7 



(112) 



as a function of p. Hence, this unified cosmological model 
has been split into a two-fluid model with energy transfer. 
The interaction Q is obtained after using the equation of 
state p = p(p), the conservation equation p' = —p — p 
and p" = p + p — p' in the source equation ©. Thus, the 
interaction 



1 

A7 



IdxP + {p + p) 1 - (7c + lx) + ~r 



(113) 

becomes a function of p. 

We conclude that instead of considering interacting 
and unified models separately, it is better to see them 
as equivalent models. However, the conclusion could be 
different when the barotropic indices 7c and 7^ are not 
constants. 



In particular, for the simplest realization of the Chaply- 
gin gas model we have investigated in the previous section 
(with 7 C = 70 and j x = 0), we obtain 



PcCh 



rch 



B 



PxCh 



B 



(118) 
(119) 



and the interaction (|116[) reduces to On the other 

hand, by inserting the Eq. into the Eqs. (|TO|) - pTg)) 
we can express the relevant quantities characterizing the 
interacting two-fluid model in terms of the scale factor. 

For expanding universes and 70(1 + a) > 0, the total 
energy density (|57|) p a+1 —¥ B/jo and the ratio has the 
non-vanishing limit rch ^ ~lx/lc (see the Eq. (|117p ). 
Thus, the above decomposition of the Chaplygin gas 
could be used to alleviate the problem of coincidence. 
However, the simplest realization of the Chaplygin gas 
(|1 19(1 has a vanishing limit for the ratio rch ~^ and it 
is unusefulness to solve the problem of coincidence. Also, 
for a = 70 = 1, we get the original version of the Chap- 
lygin gas. 



The Chaplygin gas as an interacting two-fluid 
model 



VI. CONCLUSIONS 



For illustration, we consider the class of Chaplygin 
gases generated by the equation of state ([98]) and ex- 
press it as an interacting two-fluid model. To this end, 
we insert p = p(p), given by the Eq. (|98p . into Eqs. 
(|112[) - (|113p and find the dark matter and dark energy 
densities 



P 

A7 

. _P_ 

' A7 

and the interaction Q 



PcCh - 
PxCh 



B 



To - lx ~ —fi 
B 

7c -10 + 



(114) 
(115) 



Qch = a 



IcPcPx + JxP x 



+ (1 + a)(7o - 7c )( 7c/ o c + jxPx) + IxPx 



(116) 



It can be written as a function of the total energy den- 
sity p by replacing the energy densities f|114|> - (|115[) in 
()116|) or as a function of the scale factor by using Eqs. 
(|9"T|) and (|1 14fi - (|1 16|1 . The interacting models, associated 
with the energy densities (|114[) - (|115j) and the interaction 
(|116| . are parametrized by 7 C and "f x . It means that we 
have a set of coupled two-fluid models which are related 
with a given unified model and they are described by the 
effective equation of state (|98|) . However, these models 
differ by the ratio r C h = PcCh/PxCh 



a+l 



rch 



-B + (70 - 7 x )p' 
B + (7c - 7o)/0 Q+1 



(117) 



We have investigated several models of dark matter 
and dark energy with energy transfer and shown that 
they can be considered as unified models, where both 
dark components are replaced by an effective one-fluid 
description with an effective equation of state. The two 
coupled equations describing the interacting model have 
been combined to obtain a second order differential equa- 
tion for the total energy density, "the source equation". 
We have assumed a separable interaction Q = ^(5(7,7'), 
which includes a large set of cases investigated in the 
literature, and found the conditions of stability for the 
scaling solutions. 

We have presented the "linear interaction I" and im- 
posed the conditions of stability to restrict the fourth 
constants of Qi = bip + b2p c + b3p x +b4p' . Then, we have 
found the two stationary solutions 7 S and 7/" of the evo- 
lution equation for the effective barotropic index, being 
7 S the attractor solution. Interestingly, the existence of 
the attractor solution is linked to the fulfillment of the re- 
quirement Qi(js) < 0, indicating that the energy is being 
transfcred from the dark energy to the dark matter. We 
have considered several particular examples by analyzing 
each term of Qi separately and found that some of these 
terms i.e., the total energy density, its derivative and the 
dark energy density are satisfactory couplings because in 
each case, the ratio r s = p cs j p xs is and attractor. These 
simple models may alleviate the problem of coincidence. 
Although, a coupling proportional to the dark matter en- 
ergy density does not lead to stable solutions however, it 
can work when it is combined with the remaining terms 
of Qi. 
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Taking into account that p c , p x , p' c , p' x , p, p' and p" 
can be written as linear functions of the basis elements 
p, p' and p", we have investigated the "linear interaction 
II" , by introducing a linear combination of the former 
terms Q L = Cip c + c 2 p x + c 3 p' c + c±p' x + c 5 p + c 6 p' + c 7 p" , 
and reduced it to a linear combination of the latter ones 
Ql = c'ip + c' 2 p' + c'sp" . As far as we know, this interac- 
tion has not been investigated in the literature. We have 
imposed the conditions of stability and obtained the sta- 
tionary solution for the effective barotropic index, 7 S and 
7^ , being 7 S the attractor solution, the exact scale fac- 
tor and the effective barotropic index in implicit form. 
To generalize the above coupling, we have gone a step 
more by adding a constant and a well-behaved function 
of the scale factor to the " linear interaction II" . Then, we 
have introduced the "general linear interaction", Q q l = 
7s7^A e jy/A7 + Ql, and a "more general linear interac- 
tion" Q mgL =it(ll - u)A e// / 7s A 7 + /(r?)/A 7 + Q L . 
Since the relative weight between the two first terms in 
QmgL determines the final behavior of the scale factor, 
the models generated by those interactions have given 
several alternatives to the ACDM model. 

We have presented a class on "nonlinear interaction" 
Q nL = (c'{p~ 1 p' 2 +4p+4p'+c1p"+f(r 1 )p c "/^~^)/A 1 , 
which leads to several important results. Although the 
source equation becomes a nonlinear differential equa- 
tion, we have linearized and reduced it to the equation 
of motion for a forced linear oscillator with dissipative or 
antidissipative effects. Then, the analysis of this interac- 
tion has been separated into two main parts, the homo- 
geneous case and the inhomogeneous one. In general, we 
have found that the equation of state of the effective one- 
fluid model depends explicitly on the scale factor. On the 
other hand, the universe evolves to a power-law scenario 
for large cosmological times. However for interactions 
having the form Q u l = o>"f c p c p x /p with a a constant, we 
have shown that the effective equation of state becomes 
that of the Chaplygin gas when the dark energy compo- 



nent is described by some kind of vacuum energy density 
i.e., j x = 0. Also, we have investigated the "reduced 
unified model" obtained for the particular nonlinear in- 
teraction leading to the source equation x" = 0. In this 
case we have shown that the effective equation of state 
can be interpreted as a generalization of the polytropic 
equation of state. 

Generically, when there are no restrictions on Q n L, 
the equation of state of the effective one-fluid model de- 
fines what we have called the "relaxed Chaplygin gas 
model" . It contains various generalizations of the Chap- 
lygin gas, including the variable modified Chaplygin gas 
model with equation of state p = Ap+B(a) / p a . To learn 
more about this nonlinear interaction, we have included 
a short analysis devoted to the structural stability of the 
source equation solutions. This has been made by es- 
tablishing some connection between those solutions and 
the respective function /(a) contained in the nonlinear 
interaction. For a selected set of functions /(a), we have 
shown that the final asymptotic solutions of the source 
equation can be easily obtained. 

Finally, we have given a prescription to obtain an inter- 
acting model starting form a unified one and found the 
energy densities of the dark components together with 
the interaction term that generates the former model. 
For illustration we have applied the prescription to the 
Chaplygin gas getting the relevant quantities which de- 
fine the interacting model. 
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